Phase diagram of the frustrated Hubbard model 
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The Mott-Hubbard metal-insulator transition in the paramagnetic phase of the one-band Hubbard 
model has long been used to describe similar features in real materials like V2O3. Here we show 
that this transition is hidden inside a rather robust antiferromagnetic insulator even in the presence 
of comparatively strong magnetic frustration. This result raises the question of the relevance of the 
Mott-Hubbard metal-insulator transition for the generic phase diagram of the one-band Hubbard 
model. 



The microscopic description of magnetism and metal- 
insulator transitions constitutes one of the major research 
activities in modern solid state theory. Especially tran- 
sition metal compounds like V2O3, LaTiOa, NiS2-a;Se2: 
and the cuprates show metal-insulator transitions and 
magnetic order depending on composition, pressure or 
other control parameters The paramagnetic insu- 
lating phase observed in these materials is believed to 
be a so-called Mott-Hubbard insulator due to electron- 
electron correlations; in contrast to Slater or band insu- 
lators hke SrTiOs. 

The simplest model showing both magnetism and a 
correlation-induced metal-insulator transition (MIT) is 
the one-band Hubbard model 
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Considerable progress in understanding the physics of 
this simple but nevertheless non-trivial model has been 
achieved in the last decade through the development of 
the dynamical mean- field theory (DMFT) HHli. In 
particular, the phase diagram for the unfrustrated Hub- 
bard model is very well understood 4, ^ . At half- filling 
the physics is dominated by an antiferromagnetic in- 
sulating phase (AFI) for all C/ > with a maximum 
T/v ~ 0.15W^ around U W, where W is the bandwidth 
of the non-interacting system. For finite doping, the an- 
tiferromagnetic phase persists up to a critical doping dc 
and in addition shows phase separation J^, _8]. For 
very large values of U the antiferromagnetic phase is re- 
placed by a small region of Nagaoka type ferromagnetism 

Under rather special assumptions it is possible 
to introduce complete magnetic frustration and suppress 
the antiferromagnetic phase. In this case, a transition 
from a paramagnetic metal (PM) to a paramagnetic in- 
sulator (PI) becomes visible at half filling. At T = it 
occurs at a value of the Coulomb parameter Uc ~ 1.5W 

tlyJlSj. Interestingly, the transition is of first order 
for T > with a second order end point at a 
Te « 0.017VF and Uc « 1.2W. Note that < TT^. 

A closer look at the phase diagram of V2O3 [l5| re- 
veals a strikingly similar scenario, and indeed the DMFT 



results for the Hubbard model have been used as a qual- 
itative explanation . For a proper description of 
this material, however, the antiferromagnetic phase be- 
low Tn ~ 160K has to be taken into account. It was 
argued and generally accepted j5j that the introduction of 
partial magnetic frustration will lead to the anticipated 
situation, where the MIT extends beyond the antiferro- 
magnetic phase at low temperatures. The merging of 
these two transitions presents an interesting problem on 
its own, because it is commonly believed that the mag- 
netic transition should be of second order. Furthermore, 
previous results for a system with magnetic frustration 
show, in addition to the desired effect of reducing T/v, 
an extended antiferromagnetic metallic (AFM) phase at 
small U, preceeding the transition to the AFI [l^. 
This observation suggests an appealing possibility to link 
the MIT in the paramagnetic phase with a transition 
from an AFM to an AFI. 

In this paper, we present results from a calculation us- 
ing Wilson's numerical renormalization group approach 
(NRG) [1^ and exact diagonalization techniques (ED) 
jiot to solve the DMFT equations [lEi for the Hubbard 
model with magnetic frustration at half filling. We show 
that frustration of the magnetic interactions through in- 
corporation of suitable long-range hopping does not lead 
to the previously reported sequence PM^AFM<-^AFI 
with an extended region of an AFM at T = H, E3| ■ 
Instead, we observe a first-order transition PM<-^-AFI. 
Furthermore, the reduction of Tjv is too small to result 
in the qualitative phase diagram of V2O3. 

The natural choice for studying the effect of magnetic 
frustration is the simple hypercubic lattice with nearest 
and next-nearest neighbor hopping. In the DMFT, the 
lattice structure only enters via the dispersion of the non- 
interacting band states, and the corresponding k sums 
can conveniently be transformed into energy integrals us- 
ing the free density of states (DOS) A further simpli- 
fication arises if one considers lattices in the limit of large 
coordination number. Especially for the simple hypercu- 
bic lattice the DOS then becomes a Gaussian i^J^J and 
the integrals can be performed analytically 

The investigation of magnetic properties is straight- 
forward, too. In the case of the Neel state, the lattice is 
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divided into A and B sublattices which results in a matrix 
structure of the DMFT equations |^. An antiferromag- 
netic Neel order then corresponds to a finite staggered 
magnetization mg > with — rng and to^ — —ms- 
Unfortunately, the Gaussian DOS of the hypercubic lat- 
tice has no real band edges, but stretches to infinity; the 
resulting exponential tails therefore prevent a clear dis- 
tinction between metal and insulator at T = 0, as has 
been observed in a Hartree calculation for the hyp ercu- 
bic lattice with next-nearest neighbor hopping [2l|. 

Georges et al. suggested an extension of the DMFT 
equations for the Bethe lattice which introduces magnetic 
frustration with the DOS having finite support; this re- 
sults in an analytically tractable form of the DMFT equa- 
tions even for the AB-lattice. 

For the standard Bethe lattice with nearest-neighbor 
hopping t and in the limit of infinite coordination number 
the DMFT equations on an AB-lattice acquire the form 
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Frustration is then introduced via additional terms in the 
denominators of Q (StlTrf 
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and the constraint — t\+t\. In the paramagnetic case, 
the equations © reduce to those of a standard Bethe lat- 
tice which, for example, result in the well-studied Mott- 
transition. Furthermore, despite the frustration intro- 
duced, the system is still particle-hole symmetric. Espe- 
cially for half filling this feature reduces the numerical 
effort quite drastically. Note that a similar suggestion 
exists for the hypercubic lattice 0|, too. 

Invoking the symmetry GAa{z) = GBa{z) valid for 
the Neel state, eqs. © reduce to two coupled nonlinear 
equations which we solve iteratively. In the course of the 
iterations, the quantity T,a{z) has to be calculated from 
the solution of a generalized single impurity Anderson 
model 0. For that task we employ the NRG 01, suit- 
ably extended to treat spin-polarized situations . 

Let us first discuss the results for the magnetization 
as function of U for T = 0. In the following, we fix 
t2/ti = 1/v^ « 0.58 jH and use the bandwidth W 
of the non-interacting system as our energy scale. The 
NRG results in Fig. (circles) show a completely dif- 
ferent behavior as compared to the ED data (diamonds) 
from ref..l7t Instead of a continuous increase of the stag- 



1 

0.8 - 

0.6 - 

0.4 - 

0.2 - 




r 



▲ ■ ■ A ED (6 sites) 
4 ^ Chitra et al. 

H h Hartree solution 



0.5 



U/W 



FIG. 1: Staggered magnetization ms as function of U at T = 
0. The circles are the results from NRG; the triangles from an 
ED calculation for 6 sites, while the diamonds are taken from 
ref. 11^ . For comparison, the results of a Hartree calculation 
are given by the crosses. 



gered magnetization ms for U > Uc ~ OAW as sug- 
gested by both a Hartree calculation (crosses) and the 
data from ref. llTl we find a jump in ms at a considerably 
larger Uc ~ 0.914^. To clarify this discrepancy, we per- 
formed ED calculations, resulting in the triangles shown 
in Fig.^ We find good agreement with our NRG results, 
the transition systematically approaching the NRG curve 
with increasing size of the system diagonalized in the ED 
procedure. We furthermore observed a rather strong de- 
pendence of the ED results on details of the numerical 
procedure, especially on the energy cutoff introduced in 
calculating G{z). Decreasing this cutoff systematically 
shifts the ED result towards the one found in 17]. The 
NRG, on the other hand, is stable with respect to changes 
in the parameters controlling its numerical accuracy. 
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FIG. 2: Density of states for spin up on an A lattice site as 
function of frequency. 

Another important question is the existence of an anti- 
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ferromagnetic metallic solution of the DMFT equations. 
Figure H shows the NRG results for the DOS for T = 
and spin up on an A lattice site. Due to particle-hole 
symmetry the DOS for spin down on A sites (or spin up 
on B sites) can be obtained by lu —lu. The full and 
dashed lines represent the AFI solution for U \ Uc and 
the PM solution for U Uc, respectively. Clearly, the 
magnetic solution is insulating with a well-developed gap 
at the Fermi energy. Quite generally, we were not able 
to find a stable AFM solution at T = 0. 

The discontinuity in the staggered magnetization ms 
at the transition PM<-^AFI implies a first order transition 
and the existence of a hysteresis region. Indeed, starting 




1.2 
1 

0.8 
0.6 
0.4 
0.2 



r- 


1 






i 

\ 


«m(U) from AFI 
^♦A(0) from PM 
- ^ A(0) from AFI 
«m(U)fromPM 










0.85 



0.9 



0.95 



U/W 



FIG. 3: Staggered magnetization (solid lines) and total DOS 
at the Fermi energy (dashed lines) as function of U in the 
vicinity of Uc for T = 0.0155VF. The arrows indicate that the 
DMFT solutions have been obtained by either increasing U 
(^) or decreasing U (^). 

from the paramagnet a,t U <ti Uc and increasing U results 
in a magnetization curve different from the one obtained 
by starting at U ^ Uc and decreasing U. This is appar- 
ent from Fig.|31(full lines) where a region of hysteresis can 
be observed in the staggered magnetization (for temper- 
ature T = 0.0155VF). At the same time the total DOS at 
the Fermi energy A(0) = A-^ (0) -I- (0) shows hystere- 
sis between metallic and insulating behavior in exactly 
the same U region. Note, that due to the finite temper- 
ature the DOS at the Fermi level is not exactly zero in 
the Neel state, but strongly reduced as compared to the 
metal ^V^- 

It is of course important to verify that the hysteresis 
found for small U is not some kind of artefact. This can 
most conveniently be shown by looking at the transition 
at large U. Due to the mapping of the Hubbard model 
to a Heisenberg model in this regime, one should expect 
the transition to be of second order, with the staggered 
magnetization vanishing continously like ms oc \/Tn — T 
when approaching Tjq from below. That this is indeed 
the case is apparent from Fig. ^ where we show the 



FIG. 4: Squared staggered magnetization ms as function of T 
at U/W — 2. Note that m| vanishes continously like Tjv — T 
asT\TN. 



squared staggered magnetization as function of T for 
U/W = 2. The transition is thus of second order with 
the expected mean-field exponent in this region of the 
phase diagram. 

Collecting the results for the transitions and the hys- 
teresis region for different temperatures leads to the 
phase diagram in Fig. |S1 An enlarged view of the re- 
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FIG. 5: Magnetic phase diagram for the Hubbard model with 
frustration as defined by eqs. O and t2/ti = l/\/3. The 
dotted lines inside the AFI denote the coexistence region for 
the paramagnetic MIT. The inset shows an enlarged view of 
the region with coexistence of PM and AFI. 

gion showing coexistence of PM and AFI is given in 
the inset, where the full line represents the transition 
PM^AFI with increasing U and the dashed line the 
transition AFI— s-PM with decreasing U. These two lines 
seem to merge at a value of [/ ~ IF for this particular 
value of t2 , with a critical temperature for this endpoint 
Tc « 0.02iy. Note that, even in the presence of such a 
sizeable t2, the antiferromagnetic phase still completely 
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encompasses the paramagnetic MIT (dotted lines in the 
main panel of Fig.lSlfl^'). 

It is, of course, interesting to see how the magnetic 
phase evolves with increasing t2 and in particular how 
its boundary crosses the paramagnetic MIT. We find that 
increasing t2 does not change the form of the magnetic 
phase in Fig. qualitatively, but mainly shifts the lower 
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FIG. 6: Schematic evolution of the magnetic phase diagram 
with increasing frustration. The dots on the phase transition 
lines denote the critical endpoints of the first order transitions. 

critical U and decreases the maximum Tjy. The calcu- 
lated estimates for those two quantities as function of t2 
lead to the schematic evolution of the phase diagram pre- 
sented in Fig.|SjL-c. Here, only the true phase boundaries 
are shown. A direct calculation of the free energy at fi- 
nite temperatures is presently not possible with the NRG 
method, so we cannot calculate the actual transition line 
separating the paramagnetic and AF phases. The tran- 
sition lines in Fig. |H^-c are therefore a guide to the eye 
only. For the Mott transition, the position of the actual 
transition line has been calculated in ref. |3 

Figure IHli shows the qualitative phase diagram corre- 
sponding to Fig. |31 with the line of first order transitions 
ending in a critical point. Upon further increasing the 
value of t2, the first order transition lines from both the 
PM^AFI and the Mott transition cross (Fig. EJ)), thus 
exposing a finite region of the Mott insulator and a tran- 
sition PI<-^AFI. Finally, for even higher values of t2, the 
PM^AFI transition at T = approaches the Mott tran- 
sition and T/v is reduced significantly (Fig.lHt). Note that 
in the limiting case t2 ~ ti the AFI phase completely 
vanishes due to the structure of the DMFT equations 
However, for t2 — *■ ii there is always a finite antifer- 
romagnetic exchange J oc {ti ^ ^2) /U which is sufficient 
to stabilize an antiferromagnetic ground state for U > Uc 
of the Mott transition. 

From these results we conclude that frustration as in- 
troduced via eqs. l|2Jl is not sufficient to qualitatively re- 
produce the phase diagram of materials like V2O3. In 



particular, the Mott transition extends beyond the AFI 
region only for unphysically large values of t2. 

The question remains whether it is possible at all to 
reproduce qualitatively the scenario observed in V2O3 
within a one-band model. Based on our results reported 
here, we rather believe that one has to take into ac- 
count additional degrees of freedom, for example phonons 
(within a Holstein-Hubbard model) or orbital degenera- 
cies (within a multi-band Hubbard model). 
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